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ABSTRACT 

In this Letter, the radial velocity variability of the chemically peculiar star e Ursae Ma- 
joris (e UMa) from the sharp cores of the hydrogen lines is investigated. This study is 
based on the ELODIE archival data obtained at different phases of the rotational cycle. 
The star exhibits low- amplitude radial velocity variations with a period of P = 5.0887 
d. The best Kcplerian solution yields an eccentricity e — 0.503 and a minimum mass 
'^14.7Afjup on the hypothesis that the rotational axis of e UMa is perpendicular to 
the orbital plane. This result indicate that the companion is the brown-dwarf with the 
projected semi-amplitude variation of the radial velocity K2 = 135.9 km s^^ and the 
sine of inclination times semi-major axis a2sin(i) = 0.055 au. 
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1 INTRODUCTION 

Epsilon Ursae Majoris (e UMa, HD 112185, HR 4905) is 
the brightest {V = 1.77 mag) chemically peculiar (CP) star 
and has been e xtensively studied during the last century. 
iGuthnickl l|l934l 'l established a period of 5.'^0887 d from vari- 
ations in the intensity of the C a II K line and also noti ced a 
periodic splitting of some lines. IStruve fc Hiltneill|l943t ) sub- 
sequently reported doubling of lines of Cr II, Fe II, V II and 
other elements at certain phases. Since the overall widths 
of the double lines are the same and not all lines double, 
they ruled out the orbital motion and instead suggest ed that 
the p henomenon is related to rotation of the star. iProvinI 
('1953) measured a double wave light variation with the same 
5.'*0887 period. The star is brightest when the Ca II K line 
intensity is near its minimum and other elements are near 
their maximum strength. 

Perhaps the most interesting and contr oversial aspect of 
e UM a is its the radial velocity variations. lAbt fc SnowdenI 
l|l973t) carried out the radial velocity analysis of the 62 
brightest northern CP stars for spectroscopic binaries. They 
were used special observing and measuring techniques of the 
hydrogen lines and did not find the radial velocity varia- 
tions for e UMa. Not surprisingly, they could not to measure 
the radial velocity variations because the mean internal er- 
ror per spectrum of 1.2 km s~^ is near the limit at which 
the photographic technique is useful. On the other hand. 
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[Morgan et al.l (|l978l ) detected the duplicity for the system 
e UMa with a separation of 0.053 arcsec, using speckle inter- 
ferometry technique. Although, the authors noted that the 
observed separation corresponds to the diffraction limit of 
the 2.5-m Isaac Newton telescope. 



IWoszczvk fc Jasihskil (|l980l ) measured the radial veloc- 
ity variations of many lines in the star and found sinusoidal 
variations in the radial velocities of Fe, Cr and Ti lines with 
amplitudes of about 20 km s~^ and attributed this to the ex- 
istence of a spots of enhanced abundance of these elements. 
Several surface abundance Doppler i mages of e UMa have 
been produced (IWehla.u et al.l Il982l: iRice fc Weh lau 'l990'; 
lHatzeslll99ll : lRice et al.lll997l : lHolmgren fc RicdbOOQ l . Most 
of the published Doppl er images relate to Fe , O, Ca, or 
Cr abundance. Recently. iLueftinger et al] (|2003l ) determined 
for e UMa for the first time the abundance distributions of 
Mn, Ti, Sr, and Mg. Attempt to determine orbital motion 
using these spectral lines would be impossible. However, the 
hydrogen lines do not show significant rotational effect in 
their radial velocity. They would be preferable for the mea- 
surement of binary motion if the radial velocity of the hy- 
drogen lines could be measured with precision. 

In this Letter, we present the radial velocity measure- 
ment of the hydrogen lines in the spectrum of CP star 
e UMa. The ELODIE observations of this star and data re- 
duction are described in Sect. 2. The orbital solution derived 
for the satellite candidate are presented in Sect. 3. Section 4 
reports discussion of our results for e UMa and conclusions 
are presented in Sect. 5. 
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Table 1. Measured radial velocities of t UMa. All data are relative to the Solar system barycentre. 
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Table 2. Best Keplorian orbital solution derived for e UMa 
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Value 
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5.0887 
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e 
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K 
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0.067 


/(m) 
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2 OBSERVATIONAL MATERIAL AND DATA 
REDUCTION 

The spe ctra of e UMa were retrieved from the ELODIE 
archive (|Moultaka et al.l |2004| ). This archive contains the 
complete collection of high-resolution echell e spectra using 
the ELODIE fiber-fed echelle spectrograph (|Baranne et al.l 
[l99d) mounted on the 1.93-m telescope at the Haute- 
Provence Observatory (France). The spectra have a res- 
olution (A/AA) of about 42000. The archived signal-to- 
noise ratio was between ~200 and ~400 in the spec- 
tral region near A 5550 A. In addition, eight spec- 
tra obtained 15 and 16 June 2000 at the same tele- 
scope were retrieved from the Hypercat Fits Archive 
( http: / /leda.univ-lyonl .fr /II/ spectrophot ometry.html| . In 
Table 1, each spectrum is presented by its file name, ex- 



posure time, Julian date of the observations and its corre- 
sponding phase. Note that the phases were computed with 
the ephemeris from Table 2 (see below). 



Special technique of radial velocity measurement of the 
hydrogen lines was used and should be explained. Well 
known, that in slowly rotating late B-type stars the hydro- 
gen lines have broad wings and sha rp cores. For e UMa th e 
value of V sin i is equal 35 km s^^ (|Lueftinger et al1l2003l ). 
Unfortunately, the wings of the hydrogen lines can be af- 
fected by different spectral lines. For example, in spectrum 
of e UMa the wing of Hg line affected by many Cr lines 
l|Ziznovskv fc Zverkdll995h . On the other hand, the sharp 
cores of the hydrogen lines are not affected by spectral lines 
and the flux is formed in the upper layers of the atmosphere. 
Thus our technique was to use only the sharp cores of the hy- 
drogen lines in spectra of e UMa. The spectra were processed 
using the spectral reduction software SPE developed by S. 
Sergeev at Crimean Astrophysical Observatory (CrAO). The 
program allows detecting the variations of the centre gravity 
of the sharp cores of the hydrogen lines. After processing, 
all the spectra were corrected for the motion of the Earth 
around the Sun. For example, the demonstration of the po- 
sitional variability of the core Ha line is presented in Fig. 1. 
In Table 1, each spectrum is presented by its radial velocities 
computed from the hydrogen lines, the mean radial velocity 
and errors of the mean radial velocity. 
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Figure 1. The observed intensity profiles of the cores of Ha line 
obtained on May 21, 2000 (phase = 0.93) and on June 6, 2000 
(phase = 0.07) marked by the solid and dashed lines, respectively. 
The velocity scale is given with respect to the A = 6562.797 A. 



3 ORBITAL PARAMETERS 

Orbital elements have been determined by a non- 
linear least-squares fitting of the mean radial ve- 
locities from Table 1 using the program BINARY 
writing D.H. Gudehus from Georgia State University 
|http:/ /www.chara. gsu.edu/~gudehus/binary.litmlV The 
solution for a single-lined binary is modelled by up to six 
parameters: 

(i) P period 

(ii) T time of periastron passage 

(iii) e eccentricity 

(iv) Vo system radial velocity 

(v) uj longitude of periastron 

(vi) a sin(i) sine of inclination times semi-major axis 



The expected radial velocities are 
RV = K[cos {6 + Lu) + e cos uj] 



(1) 



where 8 is the angular position of the star measured from the 
centre of mass at a given instant. The program also calculate 
the projected semi-amplitude variation of the radial velocity: 



K 



27ra sin i 



pVT 



(2) 



though this is never used as a parameter in the solution and 
the mass function: 



M| sin^ 



(Ml + M2)2 



(3) 



The program BINARY gives the estimated standard devia- 
tions of the orbital parameters as well. The orbital solution 
of Table 2 was obtained by fitting a Keplerian orbit to the 
29 ELODIE radial velocity measurements. Note that most of 
observations were obtained between JD=2451568 (January 
2000) and JD=2451714 (June 2000) (see Table 1). 

Experience shows that the best Keplerian fit to the data 
with the fixed period P — 5.0887 d, the eccentricity e = 0.503 
and the semi- amplitude K — 0.634 km s~^. The parame- 
ters of the best Keplerian orbital solution for e UMa are 
presented in Table 2. In close binary system with Bp-Ap 
stars, there is evidence for a tendency toward synchroniza- 
tion between the rotational and orbital moti ons. This effect 
is thought to be produced by the tidal forces (|Gerbaldi et al.l 



E -9.0 




-'--km 



Phase 

Figure 2. Top : Phase diagram of the radial velocity measure- 
ments and Keplerian orbital solution for e UMa. Bottom : Resid- 
uals around the solution. 



Il985h . The radial velocity curve is displayed in Fig. 2 with 
the residuals around solution. A linear trend is not observed 
in the residuals around the orbital solution that can be ex- 
plained by the absent of a second companion in a longer- 
period orbit. Although, the weighted r.m.s. around the best 
Keplerian solution (cr(O-C)) is equal to 0.131 km s~^. This 
value is a bit large compared to the typical radial velocity 
measurements fro m the ELODIE cross-correlation function 
iNaef et al.ll2003l ). 



4 DISCUSSION 

The presence of the spots on the stellar surface of e UMa 
can change the observed spectral line profiles and induces 
a periodic radial velocity signal similar to the one expected 
from the presence of a satellite. The hydrogen lines analysis 
is one of the best tools to discriminate between radial ve- 
locity variations due to changes in the spectral line shapes 
and variations due to the real orbital motion of the star. It is 
obviously of interest to compare the phase diagram of the ra- 
dial velocity of e UMa derived above with the phase diagram 
of the radial velocity computed from metallic lines. In this 
way, we selected the spectral region of the very prominent 
unblended Cr II line at A 4558 A. 

Figure 3 shows the strong photosphere lines in the spec- 
tral region of the Cr II line at A 4558 A obtained at phases 
before and after the epoch of periastron passage. The radial 
velocity shift of the Cr II AA 4554, 4558 and 4565 AA lines 
and the Fe II A 4555 A line at different phases is clearly seen. 
The Ti II line at A 4563 A shows the line doubling that ap- 
pears at phase 0.84 and disappears at phase 0.19. The max- 
imum splitting of this line is at phase 0.04. The splitting 
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Figure 3. Positional variations of strong photosphere lines in the 
spectral region of the Cr II line at A 4558 A obtained on June 
15, 2000 (phase = 0.84) and on January 26, 2000 (phase = 0.19) 
marked by the solid and dashed lines, respectively. The velocity 
scale is given with respect to the A = 4558.65 A. 



of the lines was first observed by IStruve fc Hiltnej ()l943l ) 
in the spectrum of e UMa. However, the authors concluded 
that the doubling is not caused by orbital motion but may 
be due to a combination of the physical effects with Doppler 
effe ct in a rotatin g star. 

IWadd (Il997l ) estimated the masses for 10 magnetic CP 
stars using the position of stars in the logJ?Q - logToff plane. 
Position of e UMa in the Io^Rq - logTcff plane seems to 
be among the most evolved CP stars an d gives the value 
of M eUMa = 3.O±O.4M0 . According to iLueftinger et al.l 
(2003) the radius of e UMa is equal to 4.2±0.27?q cor- 
responding to their choice of vsvai — 35 km s^^ and 
i = 45° using the trigonometric parallax measured by Hip- 
parcos, TT = 40.30 mas (ESO, 1997) and an angular diame- 
ter of 1.561 mas. The effe ctive temperatu re of e UMa was 
taken from the paper by ISokolovl l| 19981 ) and is equal to 
9340±530 K. On the hypothesis that the rotational axis of 
e UMa is perpendicular to the orbital plane we can estimate 
the mass of the secondary star. For the value of Mi = S.OMq 
and i = 45° Eq.(3) gives M2 = 0.014Mq. This resuU gives 
a value ~14.7Mjup, strongly suggesting that the companion 
is in the typical brown-dwarf regime. 

If we know the value Mi / M2 then it is possible to es- 
timate the projected semi-amplitude variation of the radial 
velocity for the companion according to formula: 



(4) 



where Ki is the projected semi-amplitude variation of the 
radial velocity for e UMa taken from Table 2. For the value 
of Ki = 0.634 km s~^ Eq.(4) gives K2 = 135.9 km s'V 
Thus, we can estimate the sine of inclination times semi- 
major axis using Eq.(2). Computation gives the value of 
a2sin(i) — 0.055 an. These estimates show that the proposed 
brown-dwarf is quite close to the surface of e UMa at peri- 
astron. But, such close orbits are not new. For example, the 
subgiant star HD 118203 have the planet with eccentric or- 
bit (e = 0.31), the period of P = 6.1335 d, and i s close to 
its parent star (a — 0.06 au) (|Da Silva et al.ll2005l ). 

Another way to interpret t he radial velocity va riations is 
the radial pulsation of e UMa. iRetter et al.l (|2004| ) are anal- 
ysed observations of e UMa obtained with the star tracker 
on the Wide Field Infrared Explorer satellite. The authors 



observed that a light curve has about 2 per ce nt ampli t ude o f 
photometric variability. On the other hand, iMolnarl (Il975l ) 
has presented ultraviolet light curves for e UMa from the 
OAO-2 satellite which indicate that the photometric varia- 
tions of this star are due to variable ultraviolet absorption 
effects which redistribute flux into the visible region (see h is 
Fig. 5). Note that no colour index changes (|ProviDl|l953l ). 
Certainly the radial pulsation appears unlikely given that 
the rotational period is synchronized with the orbital pe- 
riod of e UMa. 



5 CONCLUSIONS 

The archival ELODIE high-resolution echelle spectra of 
e UMa permit us to analyse the radial velocity variations 
of the sharp cores of the hydrogen lines. This allowed de- 
termining the orbital elements of binary system for the CP 
star e UMa. The best Keplerian fit to the data shown that 
the rotational period is synchronized with the orbital period. 
We are estimated the mass of the secondary star which is 
equal ~14.7Mjup. This result indicate that the companion 
is the brown-dwarf with the projected semi-amplitude vari- 
ation of the radial velocity K2 = 135.9 km s~^ and the sine 
of inclination times semi-major axis a2sin(i) = 0.055 au. 
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